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Abstract 

Using the representation theory of Ql(N, R) , we express the wave function of the 
GL(N,M) Toda chain, which two of us recently obtained by the Quantum Inverse 
Scattering Method, in terms of multiple integrals. The main tool is our generalization 
of the Gelfand-Zetlin method to the case of infinite-dimensional representations of 
Qt(N, R) . The interpretation of this generalized construction in terms of the coadjoint 
orbits is given and the connection with the Yangian Y(gl(N)) is discussed. We also give 
the hyperbolic Sutherland model eigenfunctions expressed in terms of integrals in the 
Gelfand-Zetlin representation. Using the example of the open Toda chain, we discuss 
the connection between the Quantum Inverse Scattering Method and Representation 
Theory. 
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1 Introduction 

In the late seventies two basic approaches to quantum integrable systems were formulated. 
Olshanetsky and Perelomov discovered a direct relation between the representation theory of 
non-compact semi-simple Lie groups and Sutherland-Calogero-Moser (SCM) systems 
The methods of Representation Theory provide a complete solution to the quantum SCM 
system. Soon after, Kostant and Kazhdan realized that more general classes of integrable 
systems, including open Toda chains, can be solved by the same method jlj, [5]. The wave 
functions of the SCM systems and of Toda chains appear as zonal spherical functions [U], 
[7j, and Whittaker functions [E]-PI], respectively. There is a generalization of the above 
method to a still more general class of integrable systems, formally corresponding to affine 
Lie algebras and affine quantum groups (see ^2], ED an d references therein). However, 
up to now this approach has not produced explicit representations of the wave functions in 
terms of integrals because of analytic difficulties in the corresponding representation theory. 

At the same time, another approach to quantum integrability emerged. In the works of 
Faddeev, Sklyanin, Kulish and Takhtajan ^l]-^7j 1 , the Quantum Inverse Scattering Method 
(QISM) was formulated and successively applied to numerous quantum integrable systems 
19 . A related recent development is connected with the Separation of Variables (SoV) 
method [2U], [21] ■ Nevertheless, the interrelations of QISM with Representation Theory 
have remained unclear. The main objective of this paper is to outline an approach towards 
establishing such a connection. 

1 See the summary of the early period of the development in |18|. 



2 



Our starting point was a new method of deriving the wave function for the GL(N, E) 
Toda chain given in [22] • This had no obvious relationship with the Representation Theory 
approach, but instead QISM was used to construct the wave function. The understanding of 
this construction in terms of Representation Theory, and the comparison of these two basic 
approaches to quantum integrability for this simple integrable system, should provide addi- 
tional insight into the nature of QISM. In this paper we give an interpretation of the integral 
expression of the wave function of the open Toda chain, obtained in (221; * n terms of the 
Representation Theory approach to quantum integrability. The basic tool is a generalization 
of the Gelfand-Zetlin method [23], |23) to the case of infinite-dimensional representations of 
gl(N, E) . This gives a construction of a representation of U(gl(N, E)) in terms of difference 
operators 2 , which one can use for finding explicit expressions for Whittaker vectors. Thus, 
we reproduce the integral formula for Whittaker functions obtained previously in |22j by 
QISM. 

Comparing the methods of [22J with those proposed in this paper reveals a deep connec- 
tion between the SoV [2D] and Gelfand-Zetlin type representations. In j2Ej, |2D] the separated 
variables for the classical Toda chain were found. They could be successfully quantized to 
give an explicit solution to the corresponding quantum theory. In this paper we use the 
generalization of the Gelfand-Zetlin representation for gl(N, E) , which gives rise to a much 
larger set of separated variables ( |iV(iV — 1) instead of (N — 1) ). The idea of introducing 
the enlarged set of separated variables can be traced back to the work of Gelfand and Ki- 
rillov on the field of fractions of the universal enveloping algebra [2Z] ■ On the other hand, the 
representation which we have constructed appears to be closely connected with the represen- 
tation theory of the Yangian Y(q[(N)) . The enlargement of the set of variables for solving 
the quantum problem is a crucial part of the Representation Theory approach to quantum 
integrability. One may hope that understanding the correct notion of separated variables in 
Representation Theory will turn out to be fruitful for both the theory of Quantum Integrable 
Systems and Representation Theory. 

The plan of the paper is as follows. In section 2 we construct the main tool we will be 
using throughout the paper: we give a direct generalization of the Gelfand-Zetlin method 
to the case of infinite-dimensional representations of U(gl(N)) . The crucial point is to 
use continuous unconstrained variables instead of integer-valued Gelfand-Zetlin parameters. 
Another type of analytical continuation was considered in [21] (see also |2E1)- The main 
results of this section are the analytic construction of a representation of U(gl(N)) and 
explicit description of mutually dual Whittaker modules. 

In section 3 we give the interpretation of the representation described in section 2 in terms 
of the Kirillov-Kostant orbit method [2H]. We give an explicit parameterization of the finite 
cover of the open part of the coadjoint orbit of GL(N, E) with stabilizer H N = GL(1, E)®^ . 
We show that this parameterization yields Darboux coordinates with respect to the canonical 
Kirillov-Kostant symplectic structure on the orbit. The above construction is a version of 
a construction by Alekseev, Faddeev and Shatashvili [30], [3T| adopted for non-compact 
groups. Our parameterization yields the classical counterparts for the generators of gl(N) 

incidentally, a version of this construction for the afhne algebra si(2) k has been discovered previously 
in [25]. The generalization to the case sl(N) k will be given elsewhere. 
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given in section 2. A subtle point of the orbit method interpretation is that we use only 
the finite cover of the open part of the orbit. This leads to a representation of Ql(N, K) 
which cannot be integrated to a representation of the group. However, it is obvious that the 
action of the Cartan subalgebra integrates to the action of a Cartan torus which is sufficient 
for constructing solutions of the open Toda chain. We end this section by giving contour 
integral formulas for the generators of Ql(N, R) . 

In section 4 we relate the Gelfand-Zetlin parameterization to the representation theory 
of the Yangian [32] • (Closely connected ideas for the classical Gelfand-Zetlin construction 
were developed in [33 -[35J). In this section we derive explicit formulas in the Gelfand- 
Zetlin parameterization for the images of the Drinfeld generators *4. n (A), £> n (A), C„(A) under 
the natural homomorphism: Y(qI(N)) — > U(qI(N)) . This result will be used in section 
6 to establish the connection of the Gelfand-Zetlin parameterization with the QISM-based 
approach to the integrability of the quantum open Toda chain. 

In section 5 we prove two of the main results of this paper, Theorems 5.1 and 5.2, 
by using the methods developed in section 2. We apply the representation constructed in 
section 2 to find the explicit solutions of the quantum open Toda chain and the quantum 
hyperbolic Sutherland model. For the Toda chain, we reproduce the results of (221- This 
bridges the gap between the Representation Theory and QISM approaches for these theories. 
The corresponding integral representation for the wave functions of the quantum hyperbolic 
Sutherland model appears to be new. 

In the last section, we compare the QISM-based approach to the quantum integrability 
of the open Toda chain to the Representation Theory approach. 

Acknowledgments: We would like to thank the Max Planck Institut fiir Mathematik 
(Bonn, Germany), where the paper was finished, for financial support and its stimulating 
atmosphere. 

The authors are grateful to V. Gorbunov, A. Levin, Y. Neretin, M. Olshanetsky, and M. 
Semenov-Tian-Shansky for useful discussions. 
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2 Whittaker modules in the Gelfand-Zetlin represen- 
tation 

2.1 Generalization of the Gelfand-Zetlin construction to infinite- 
dimensional representations 

We start by recalling the original Gelfand-Zetlin construction for the canonical basis in the 
space of an irreducible finite-dimensional representation of GL(N,C) [2B], 121] ■ It is based 
on the two well-known facts. 

Any irreducible finite-dimensional representation is uniquely determined by an n -dimen- 
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sional vector m N = (m N i, . . . , m NN ) with integer entries such that m N \ > . . . > m NN . The 
vector m N is called the label or the highest weight of the representation. 

Consider GL(N — 1,C) as the subgroup of GL(N,C) embedded in standard way. Then, 
each of the irreducible representations of the subgroup GL(N — 1, C) enters into the decom- 
position of a finite-dimensional irreducible representation of GL(N, C) with multiplicity not 
greater than one. If we denote by m N -i the labels of the representations involved in the 
decomposition, then their entries satisfy the additional conditions: m N j > m N -ij > m N j + i , 
for all j = l,...,N- 1. 

If we continue the decomposition step by step, we obtain a decomposition of an irre- 
ducible finite-dimensional representation of GL(N,G) into a direct sum of one-dimensional 
subspaces. Fixing an element in each one- dimensional summand yields a Gelfand-Zetlin 
basis. The vectors in this basis are labeled by triangular arrays, consisting of all m n , 
n — I, . . . , N , with entries satisfying the above mentioned constraints. 

Now, we wish to generalize the Gelfand-Zetlin (GZ) construction to infinite-dimensional 
representations of the universal enveloping algebra U(gl(N)) . 

We start with the following generalization of the Gelfand-Zetlin representation of fll(iV) . 
Namely, let ("f 1 . . . ~f N ) be a triangular array consisting of |iV(iV + 1) variables -f n = 

(inl, Inn) G C n ; U — 1, . . . , N . 

Proposition 2.1 Let M be the space of meromorphic functions in \N{N — 1) variables 
7 1; . . . , "y N _i . Then the operators 3 

- n n—l 

E nn = ( ^ Inj - E Tn-ij) ! (ri = 1, . . . , N), (2.1a) 

i=i 3=1 

n+l 

i n EI (7nj - 7n+l,r - f ) 
£ . . = __L V" 1=1 e - ind ^ . ( n = i N — l) (2-lb) 

^ ih^ minj -Ins) 

3 s^j 



n-1 

i n EI (inj ~ 7n-V + f ) 



S n+1 , n = ± J] ^TTT S ^ ! (n = 1, . . . , - 1) ( 2 - lc ) 

^ ^ 11 (7nj - 7ns) 

/orm a representation of g[(JV) m M . 

Proof. It is sufficient to show that the generators (|2.1jl satisfy the standard commutation 
relations 

[-S'rarj) -^m,TO+l] i$nm ( ^n,m+l)-£'m,m+l j 
[-^nni -Em+l,m] (^nm ( ^n,m+l)-^'m+l,nu (2-2) 

[-^n/ra+lj -^m+l,m] (-^nn -£/ n -(_i irt 4.i)5 rlm 



3 Actually, we construct the family of representations depending on the auxiliary parameter h . One could 
rescale the variables j n j — > fvy n j to get rid of this parameter. However, we will retain h in the formulas to 
make the connection with standard notations in the theory of quantum integrable systems more explicit. 
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and obey the Serre relations 




(2.3) 



These follow by direct verification. By defining the composite root generators recursively 
as Ej k = [Ej m , E mk ] for j < m < k and j > m > k , one obtains the complete set of 
commutation relations 



The representation described above is parameterized by the label *y N = (7^1, • • • ,7jvjv) 
and extends to a representation of the universal enveloping algebra U(gl(N)) . Obviously 
there is a natural representation of any subalgebra U(g[(n)) C U(gl(N)), n = 1, . . . , N — 1 
on the same space M with the label 7 n = (j n i, . . . , 7„„) . Let Z(gl(n)) be the centre of 



Definition 2.1 We say that a U(gl(n)) -module V admits an infinitesimal character £ if 
there is a homomorphism £ : Z(gl(n)) — > C such that zv = £(z)v for all z&Z(gl(n)),v G V. 

It is possible to show that the U(gi(N)) -module M defined above admits an infinitesimal 
character. Actually, a more general statement holds: 

Proposition 2.2 Each central element of U(gl(n)) acts on M via multiplication by a sym- 
metric polynomial in the variables . 

Proof. Let us use the induction over n . For n — 1 the statement is obvious. Notice that 
the operators u G U(gl(m)) acting on M are invariant under permutations in each subset of 
variables 7 1? . . . , 7 m . Now suppose that the action of the centres Z(gl(m)), m— 1, . . . , n — 1 
are given by the functionally independent symmetric polynomials of 7 m , m — 1, . . . ,n — 1. 
Let us prove that the centres of Z(Q{(n)) act on M analogously. The elements of the centre 
Z(gl(n)) commute with Z(gl(m)) and thus do not depend on the shifts over j m . On the 
other hand, the elements of Z(gl(n)) commute with all root generators E^, i,j — 1, . . . ,n 
and, hence, are ih- periodic with respect to any variable 7 mj -, m = l,...,n — 1. The 
polynomial generators of the centre of U(gl(n)) may act only as rational functions in this 
representation and thus we conclude that they do not depend on the variables 7 mj -, m = 
l,...,n — 1. Taking into account that the generators of U(gl(n)) depend polynomially 
on 7 n we conclude that the generators of the centre Z(g[(n)) act on M via symmetric 
polynomials in j n . The functional independence may be directly verified on a particular 
subset of M. ■ 



[Ejk, E lm ] — Ej rn 5ik — E lk 5. 



(2.4) 



W(fll(n)). 
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Hence, the elements z G Z($l(n)) act on the U(Ql(n)) -module M as operators of 
scalar multiplication, and the U(qI(ji)) -module M, n = 1, . . . , N , admits an infinitesimal 
character. 

In the next subsection we calculate the explicit action of the central elements on M 
using the notion of Whittaker vectors. 



2.2 Whittaker modules 



Let us now give a construction for Whittaker modules using the representation of U(qI(N)) 
described above. We first recall some facts from 

Let n + and n_ be the subalgebras of qI(N) generated, respectively, by positive and 
negative root generators. The homomorphisms (characters) x+ '■ n + ~ > Q X- '■ n - ~ > C 
are uniquely determined by their values on the simple root generators, and are called non- 
singular if the (complex) numbers x+(E nn+ i) and X-(E n+ i n ) are non-zero for all n = 
1 V- 1. 

Let V be any U = U(qI(N)) -module. Denote the action of u G U on v G V by uv . A 
vector w G V is called a Whittaker vector with respect to the character x+ if 

E n>n+1 w = x+{E n!n+ i)w , (n = 1, . . . ,N - 1), (2.5) 

and an element w' G V is called a Whittaker vector with respect to the character X- if 

E n+1 , n w' = X -(En+i,nW , (n = l,...,N-l). (2.6) 



A Whittaker vector is cyclic for V if Uw = V , and a U -module is a Whittaker module 
if it contains a cyclic Whittaker vector. The U -modules V and V are called dual if there 
exists a non-degenerate pairing (. , .) : V x V ^ C such that (Xv',v) = —(v',Xv) for all 
v G V, v' G V and X G fl[(iV) . 

We proceed with explicit formulas for Whittaker vectors corresponding to the represen- 
tation given by (j2.1|) . 



Proposition 2.3 The equations 

E n+ i in w' N = -ifT l w' N , (2.7) 



E nin+ iw N = -ih 1 w N (2.1 
for all n = 1, . . . , iV — 1 , admit the solutions 



-fEV-DE^^ 1 f (2-9) 
w N = e i- 1 [[ s n (7 n ,7 n+1 ), 

n=l 
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where 



n n+1 



Sn(7n,7n+l) =11 H & ^ 2l \ +2]' 

k=l m=l 

Proof. The following equality holds: 



n-l 



n ll (inj ~ 1n-l,r ~ f ) 

-E«,n+iWiv = -ifT l w N V" T "~ 1 , 7 • ( 2 - n ) 

^ U. VYnj - Ins) 



Using the identity 



tiy 1 . . . Jb . 



? II (•'•,''- -'V) ~ ^ -1 •••-n (2.12) 

.7=1 1 A V 5 7 fc.>0 



ki + ...+kn=m— n+1 



one arrives at (|2.8j) . Equation (|2.7|) is proved similarly. 



Evidently, the solutions (|2.9jl . (jUTUj) are not unique. The set of Whittaker vectors is closed 
under multiplication by an arbitrary ih -periodic function in the variables 7„j . Hence, there 
are infinitely many invariant subspaces in M corresponding to infinitely many Whittaker 
vectors. 

To construct irreducible submodules, let us introduce the Whittaker modules W and 
W, generated cyclically by the Whittaker vectors w N and w' N , respectively. 

Theorem 2.1 Let m n = (m n i, . . . , m nn ) be the set of non-negative integers. The Whittaker 
module W = Uw N is spanned by the elements 

N-l n 

Wm 1 ,...,m JV _ 1 = n Yl^T n "hn)w N , (2.13) 
n=l fc=l 

where o-k{~Y n ) is the elementary symmetric function of the variables j n i, . . . , 7 nn of order k: 

°k{l n )= Yl 7nh---lni k - (2.14) 

h<—<ik 

Similarly, the Whittaker module W' = Uw' N is spanned by the polynomials 

N-l n 

«'m 1 ,..,m w _ 1 = II 11^(7 J- (2-15) 

n=l k=l 

The Whittaker modules W and W' are irreducible. 
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Proof. Let us prove the statement for the module W . Using the identity e 1 ^ 1 ^ <Jk{l n ) — 
fc-i 

<Jk(-y n ) +ihYl (~^) r,T k-r-i(ln)lnj > an d the explicit formula for the action of -E n +i,n on the 

r=l 

Whittaker vector: 

71+1 

n EI ilnj ~ 7n+l,r + f ) 
^ n +l,n% = ift~ l W N ^ ^ ff7 \ , ( 2 - 16 ) 

J 

one shows, due to ()2.12|h that the Whittaker module is spanned by elements of the form 
(EH- 

To prove the irreducibility of Whittaker modules we need the following fact (0, Theorem 
3.6.2). Let M be any U -module which admits an infinitesimal character. Assume w G M 
is a Whittaker vector. Then the submodule Uw C M is irreducible. 

Proposition 12.21 and the above theorem implies that the modules W and W are irre- 
ducible. ■ 



Let us note that for any subalgebra U(gl(n)) C U(gi(N)) , 2 < n < N , the module over 

71—1 S 

the ring of the polynomials in ~f n with the basis fl II a Z sk ^Ci s) w n is a U(gl(n)) Whittaker 

s=l k=l 

module. 

Now one can calculate the explicit form of the action of the central elements of U(gl(n)) 
on the space M . It is well known that the generating function A n (X) of the central elements 
of U(gl(n)) (the Casimir operators) can be represented as follows [3*%] : 



A n {\) 

signp (A - )5 p (i),i - ... (A - ihp^Sp^n ~ ihE^ n 

P&Pn 



(2.17) 



where = |(n — 2k + 1), = l,...,n and the summation is over elements of the 
permutation group P n . 

Proposition 2.4 The operators \2.11 ) have the following form on M : 



Ai(A) = JJ(A - 7 nj ), (n = l,...,A0- 



(2.18) 



Proof. By Proposition 12. 21 it is sufficient to calculate the action of *A n (A) on the Whittaker 
vector w N (or w' N ). Due to (|2.8|1 . one has 



Ai(A)Wjv = - f )[A - - zft£„„K 

ft 

— ihy^An-kjX — l -Y 1 )E n , n _ k+ iw N . 



(2.19) 



k=2 



Furthermore, using the relation (j2.16|) it is easy to calculate the action of the composite 
generators: 



n— 1 n—k+1 k 



-1 EI (ln~p.jp 1n~p+l,r p + 2 ) 



E n , n _ k+1 w N = <r VE- E II rp ^" n u s — • ( 2 - 20 ) 

j 1= l j fc =i p =i 11 Im-W'j. Jn-p,sp) 

By substituting (j2.20j) into (|2.19|) one arrives, after some combinatoric, at ()2.18jl . Obviously, 
we can also do the analogous calculation for the Whittaker module W . ■ 

It remains to construct a pairing between W and W , and to prove that the Whittaker 
modules W and W are dual with respect to this pairing. 

Definition 2.2 Let <fi eW and ip G W . Define the pairing (,) : W <S)W ^ C by 

r N-l 

{M)= j ^0(7)0(7)^(7)]^^' (2.21) 

N(N-l) r >= 1 
R Z 3<n 

where 

N-l 

Mo(7) = II 11^ - 7n P )(e^ - e 2 ^). (2.22) 

n=2 s<p 



The integral (j2.21|) converges absolutely. Actually, this is a simple corollary of a more 
general statement. Denote the integrand in the right hand side of (|2.21j) by I N (7) . 

Lemma 2.1 For 7 „j , suc/i that min{Im7 n ,} > max{Im7„ + i fc} — I , the following inequality 

j k 

holds: 

|/zr(7)l < 1^(7)1 exp { - £ £ |Re 7nj -|} (2.23) 

^ n=l j'=l 

/or some polynomial P(i) in r y 1 , . . . , 7 JV _ 1 . 

Proof. For any x > the inequality |r(x + z?/)| < r(x)|p :r .(?/)| cosh -5 (7^) is valid for some 
polynomial p x (y) with degree depending on x . Hence, the integrand satisfies: 

|J*(7)I < \P{l)\e^\ (2.24) 



where 



N-l 



n=2 s<p 

(2.25) 

l Re 7nj - Re7 n+ i jfc | - — |Re7 JV _i 

n=l j,k j=l 
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Using the inequality 

m m m 

\ a i ~ bk \~Yl \ a i ~ a k\-^2\ b j- b k\>0, (2.26) 

j,k=l j<k j<k 

which holds for any set of real parameters {a,j,bj}; j = 1, . . . ,m (see (11.118) in [HH]), one 
proves the recursive relation S N < 2J ^_ 3 (ffjv-i — J2j |R- e 7jv-i,jl) ■ This implies ()2.23|) . ■ 



Proposition 2.5 Let ~y N G WL N . The Whittaker modules W and W are dual with respect 
to the pairing defined by \2.21\) . I.e. for any G W and if) G W, the generators X G 
qI(N, M.) possess the property 

(0,X^> = - (2-27) 

Proof. For example, consider (0, E ntn+ i%p) , where E n>n+ i is defined by (pThjl . It is easy 
to see that the expression _E , n;n+1 0-?/> does not have poles in the upper hyper-plane and it 
is therefore possible to deform the integration contour. Consider the shifts ■jnj — > + iH 
and note the difference equation 

e ^/i (7) = /i (7) [[ . (2.28) 

S9 Lj Inj Ins 

Then, using a deformation of the contour and the estimates (|2.23j) . one obtains (J2.27|) . ■ 



3 The orbit method for GL(N, R) and the Gelfand- 
Zetlin representation 

In this section we discuss the orbit method interpretation of the explicit algebraic construc- 
tion of the infinite-dimensional representation of U(gi(N)) given in the previous section. 
Similar results for the case of a compact group were obtained in jSH], [3T] . 

Let G be a Lie group, q the corresponding Lie algebra, and let q* be its dual. The 
coadjoint orbits O of G are equipped with a canonical symplectic structure given by the 
Kirillov-Kostant two-form [221 • The space of linear functions on g* is closed under the 
Kirillov-Kostant Poisson bracket and coincides with q as a Lie algebra. Consider the coad- 
joint orbits of GL(N,M.) with stabilizer H N = GL(1,M.)® N . Below, we construct Darboux 
coordinates on the finite cover of the open part of such an orbit O . We also 

give expressions in these coordinates for the restriction onto of the linear functions 

on q* . The corresponding Hamiltonian vector fields act on . The explicit realization 
of U(gl(N)) by difference operators given in the previous section may be considered as a 
"quantization" of this Poisson algebra. Let us stress that this construction gives a repre- 
sentation of U(qI(N)) which cannot be integrated to a representation of the group. Our 
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discussion below shows that this is quite natural because part of the orbit is not stable 
with respect to the action of the group on the orbit O . 

Let G = GL(N, R) and g = gl(N, R) . We identify g and its dual g* via the Killing 
form. Consider the coadjoint orbit O of the diagonal element u N G Q* with pairwise 
different diagonal entries: 7^ = (7^1, 7^2, . . . , ^ NN ) : 



(iNl ■■■ \ 



u f 



\ • • • 7JVJV / 



(3.1) 



The stabilizer of u N is H N = GL(1, R)® N and the general element of the orbit may be 
parameterized as u = g~ x u N g , where g G GL(N, M) is defined up to the left action of H N . 
There is a canonical symplectic two-form on O , given in this parameterization by 



Q = Tru N (dgg r ) 2 . 
Consider the subspace in GL(N,R) which consists of the elements: 

9 = 9n9n-i •••92. 



(3.2) 



(3.3) 



The matrices g n differ from the unit matrix only in the upper-left n x n corner and are 
defined recursively as follows. Denote the upper-left (n x n) sub-matrix of g n by f n . There 
should exist certain 7 n -ij , j 



f n 



Til 







1, . . . , n — 1 so that the following condition holds: 
°\ hn-1,1 ••• A 



V 



7n,n-l 
Inn) 



fn 



\ 





* 



: 

7n-l,n-l * 
* * 



(3-4) 



This means that the diagonal matrix conjugated by f n can have non-zero, non-diagonal 
entries only in the last column and on the last row. The matrix f n is defined uniquely by 
7 n j,7„_i i fc up to the left action by the diagonal matrix and the right action by the diagonal 
matrix with the last diagonal element being equal to one. Let us choose the following 
representatives for f n in the coset of G by the left action by diagonal matrices: 



(fnjjk — 



Q 



n—l,k 



n-1 

EI (inj - 7n-l,r) 
r=l 



Inj ~ 7n-l,fc II (inj ~ Ins) 



(k= l,...,n- 1), 



(fn)jn 



n-1 

II (inj 
r=l 



(3.5) 



In— l,r) 



EI {inj Ins) 



where Q n -i,k , k = l,...,n — 1 are additional coordinates on the orbit. This gives a 
map from R N ( N ~ 1 \A to the open part of the orbit. Here, A is a union of three 
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subspaces of RW-V : A = D x U D 2 U D 3 . The subspace L>i is a subspace where at least 
two 7 -coordinates of the same level coincide: 7^ = ^ n k- The subspace D 2 is defined 
as a subspace where at least two of the 7 -coordinates from consecutive levels coincide 
Inj = 7n-i,fc • Finally, the subspace D 3 is a subspace where at least one of Q n j is zero. 
Now, let us construct the inverse map. Consider the element of the orbit u = g~ l u N g shifted 
by the N x N unit matrix / multiplied by a formal variable A : 

T(\) = \I-u. (3.6) 

The n -th principal minor a„(A) of T(A) does not depend on g m , m < n+1 . For the general 
point of the orbit, the roots of the minors may be complex. Consider the open part of 
the orbit O such that all upper-left sub-matrices are diagonalizable and their determinants 
have distinct real roots. We also impose the condition that there are no identical roots of 
minors with consecutive ranks. Then, define the 7 -variables to be the roots of the minors: 

n 

a»(A) = ~ 7*) • (3-7) 

i=i 

Note that, by this condition the roots are defined only up to the action of the symmetric 
group S n . As a result we obtain the coordinates on the finite cover of an open part of the 
orbit. 

Now, we give the explicit expressions for the coordinates Q n j . Let b n (X) be the deter- 
minant of the n x n corner sub-matrix of the matrix obtained from T(A) by interchanging 
the n-th and the (n + l)-th columns. I.e. b n (X) = det | |T(A)5 , „ in+ i| | , where the N x N 

nxn ' 

matrix S^n+i is a unit matrix with the diagonal 2x2 sub-matrix in the n, n+1 columns 
and the n, n + 1 rows replaced by the matrix ^ ^ y ^^ mnar ^' we define c n (A) as the 

determinant of the nxn corner sub-matrix obtained from T(A) by interchanging the n-th 
and the (n+l)-th rows; i.e. c n (A) = det ||S , riiTl+ iT(A)|| . Then, the rest of the variables may 

nxn ' 

be introduced as follows: 

Qnj = b~ l (j nj )a n+1 (^ n j) = -Cn{lnj)an-l{lnj) , (3.8) 

where the last equality can be proved by a straightforward check. It is not difficult to verify 
that this definition is compatible with the explicit parameterization in terms of {^njiQnj) 
given above. Due to the condition (J3.4|) . the n-th left-upper corner sub-matrix of the matrix 
u n = g~+-L ■ ■ ■ g~ l u N g N ■ ■ ■ g n +i is diagonal matrix; a„(A) is given by ()3.7|) . The determinant 
6 n (A) is invariant under conjugation by the matrix which differs from the unit matrix only 
in the upper-left (n — 1) x (n — 1) corner. Thus, b n (\) may be represented as: 

fo n (A) = det^ 1 ) det ||(A - u n )g n S n , n+x \\ . (3.9) 

nxn v ' 

A straightforward calculation gives 

n n+1 , 

bn(X) = £ Qn] 11(7* - ^) II ■ ( 3 - 10 ) 

■ 11 _l ■ Inj ins 

3=1 r=l s^j 1 J 
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Similarly, for c n (A) one has 

n n—l , 

Cn(X) = ~ £ Qni II(7n, " 7»-l,r) ]I I ^ " ( 3 - U ) 
■ -11 / • 7n? 7ns 

This gives the inverse transformation from the finite cover of the open part of the orbit 

N(N-l) , 

to R~^—\A. 

A simple calculation of the Kirillov-Kostant symplectic form in the coordinates (7^, Q n j) 
gives the following: 

Proposition 3.1 1. The coordinates {^37 Qnj) parameterize the finite cover O^ > of the 
open part of the orbit . The product of the symmetric groups S = Yln=2 $n acts 
freely on the fibres of the projection: ir : — > and the factor is . 

2. The lift of the 2-form Vt on has the canonical form 

N-l n 

n * Q = Y,Yl d ^nj A Qn]dQny (3.12) 
n=l j=l 

The coordinates (7^, Q n j) give Darboux coordinates on the finite cover of the open 
part of the coadjoint orbit of GL(N,M) with the stabilizer conjugated to H N = 

GL(1,R)* N . 

Now let us give explicit expressions for the elements on the super-diagonal, diagonal and 
over-diagonal of the matrix u = g~~ 1 u N g in the GZ-coordinates. These give the classical 
counterparts of the generators (|2.1|) . 

Proposition 3.2 The linear functions u nn , w n ,n+i , Un+i,n on g[(N,M)* , being restricted to 
and lifted by the projection map n : — > to , have the following form in 
GZ-coordinates 4 : 

n n—l 

V"nn ^ ^ 7nj ^ ^ In— X,j > 
3=1 3=1 



n+1 

n 11 {inj - ln+l,r] 



^n,n+l 



ir = n(7n;-7n S ) Q ~ nh ( 3 - 13 ) 



3 * 



n-l 



n II (7nj - 7 



n— l,r y 



n+l,n /J t-t / s Qnj ■ 

1 11 (7nj - 7n.) 



4 There is another representation with "square roots" which is closer to that given in |2.'ij. It may be 
obtained by a different choice of parameterization for the solution of (|3.4|) . 
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Proof. Let us start with the diagonal elements. The coadjoint action of the diagonal matrix 
$ = diag(e < ^ 1 , • • • e^ N ) on u may be represented as the right action on the element g , (|3.3j) . 
It is easy to see that this is equivalent to a shift of the variables: 



(1 — > D .p4>n-<Pn+l 

(3.14) 

QnJ * QnjC N ■ 

Thus, the corresponding Hamiltonian functions u nn are: 

n n— 1 

Unn = 22 7 ™J ~ ^-1,] • ( 3 - 15 ) 

i=i i=i 

This result could also be derived from the explicit representation of the matrix u . Moreover, 
from the structure of the matrices u n , f n one could easily obtain expressions for the other 
generators. For example, 

n-l 

n n Yl \ fnj 7ra— l,r) 

Un+l,n ^ ^ (^n)n+l,j(/n)jn ^ ^ tT? ^ Qnj • (3.16) 



Quite similarly, 



, 11 (Tnj - 7n S ) 



n+1 

i 11 (7»i - 7n+l,r) 



u ^ 1= _V — O" 1 ( 3 - 17 ) 

£ n(7ni-7n,) 



There is a representation for w n n , and u n+ i in in terms of contour integrals using 

the polynomials a ra (A),6„(A), and c n (A) : 

a n (A) dA 



27ri / a n _i(A) A 



w n ,n+i = j> a n 1 (X)b n (X)dX , (3.18) 

u n +i,n = -— j> c n (X)a n 1 (X)dX . 

In the next section we discuss the connection of this representation with Drinfeld's "new 
realization" of the Yangian 



4 The Yangian Y($l(N)) and the Gelfand-Zetlin repre- 
sentation 

In the previous section we outlined the construction of the Gelfand-Zetlin parameterization 
of the coadjoint orbit of GL(N, R) . The polynomials a n (A), b n (X) and c„(A) , which appear 
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in the invariant formulation of this parameterization, bear an obvious similarity to the basic 
ingredients of Drinfeld's "new realization" of the Yangian. In this section, we derive the 
explicit expressions for the Drinfeld generators of the Yangian Y(q[(N)) in the Gelfand- 
Zetlin realization 5 . These explicit expressions will be important for a new interpretation of 
the QISM approach to the solution of the open Toda chain [2*2"j . 

We recall some well-known facts about the Yangian Y(qI(N)) , [321 ( see a ^ so the recent 
review jSH])- The Yangian Y(q[(N)) is an associative Hopf algebra generated by the elements 
T^j , where i,j = l,...,N and r = 0, . . . , oo , subject the following relations. Consider the 
N x N matrix T(A) = ||Tjj(A)||f 7=1 with operator valued entries 

oo 

T l] (X) = X5 l , + Y,T^X- r . (4.1) 

r=0 

Let 

R N (X) = I ® I + ihP/X ; P ikJl = 5 u S kj (4.2) 
be an iV 2 x N 2 numerical matrix (the Yang R -matrix). Then the relations between the 

(r) 

generators TA can be written in the standard form 

R N (X - p)(T(X) ® I) (I ® T(ji)) = (I ® T(ji))(T(X) ® I)R N (X - p). (4.3) 

The centre of the Yangian is generated by the coefficients of the following formal Laurent 
series (the quantum determinant of T(A) in the sense of [19]): 

det g T(A) 

= ^ si & n P T p(i),i( A ^ ihp[ N) ) ■ ■ ■ T p{k)ik (X - ihp[ N) ) . . . T p(N)iN (X - ihp^) , ( ' ' ' 

pePjv 

where pn = — 2n + 1), (n = 1, . . . , N) and the summation is over elements of the 
permutation group P N . Let X(X) = ||^ij(A)||" J=1 be an n x n sub-matrix of the matrix 
||Tij(A)||^ =1 . It is obvious from the explicit form of R N (X) that this sub-matrix satisfies 
an analogue of the relations (J4.3|) . The quantum determinant det g X(A) is defined similarly 
to (|4.4j) (with the evident change N — > n). 

The following way to describe the Yangian Y(gl(N)) was introduced in [32]. Let A n (A) , 
n — 1, . . . , N , be the quantum determinants of the sub- matrices, determined by the first n 
rows and columns, and let the operators B n (A), C n (A) , n = 1, . . . , iV — 1 , be the quantum 
determinants of the sub-matrices with elements Tij(X) , where i = 1, . . . , n ; j = 1, . . . ,n — 
1, n + 1 and i = 1, . . . ,n — 1, n + 1 ; j = l,...,n, respectively. The expansion coefficients of 
A n (A), B n (A), C n (A) , n — 1, . . . , N — 1 , with respect to A , together with those of Ajv(A) , 
generate the algebra Y(gl(N)) . The parts of the relations between the Drinfeld generators, 

5 For the connection of the finite-dimensional Gclfand-Zctlin representation with the Yangian see [6'6\ . 
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which we use below, are as follows: 

[A n (A), A m (/i)] = ; (n,m = l,...,N), 
[B n (A),B m (/x)] = 0; [C n (A),C m (^)] = 0; (m^n±l), 

(4.5) 

(A - // + ifi) A n (A)B n (//) = (A-^)B n (/x)An(A)+z/iA n (/i)B n (A), 
(A - fi + ift)A n (//)C„(A) = (A- / u)C„(A)A„(/i) + ^A„(A)C„(/i). 

Let A(g[(iV)) be the commutative subalgebra of Y(gl(N)) generated by A n (A), n = 
1, . . . , N . It was proved in |37| that A(qI(N)) is the maximal commutative subalgebra of 
Y(gl(N)). 

There is a natural epimorphism 71^ : Y(gl(N)) — > U(gl(N)) , 

Tc N {T jk (X)) = X5 jk - ifrE jk , (j,k = l,...,N). (4.6) 

Denote the images under n N of the generators A n (A),B n (A) , and C n (A) by A n (X),B n (\) 
and C n (\) , respectively. Let us describe the images of the Drinfeld generators. The images 
of the generators of A(q[(N)) under the homomorphism ()4.6|) have the form (j2.17|) . In 
particular, the image of A Ar (A) , which is a polynomial in A of order N, gives a generating 
function for the generators of the centre Z C U(qI(N)) . For the other generators B n (X) , 
C n (\) we have 

Lemma 4.1 

B n (X) = [A»(A), E n>n+1 ], (4.7) 
C n (A) = [E n+1<n ,A n (X)], (4.8) 

where n = 1, . . . , N — 1 . 

Proof. A direct computation using f)4.4j) and the explicit expressions for the images of the 
generators A n (A),B n (A) and C„,(A) under the homomorphism Tt N . ■ 

Using ()4.7j) . ()4.8|) and the explicit expressions for the generators Ej k introduced in Propo- 
sition 2.1 we arrive at the following 

Proposition 4.1 The operators 

n 

A(A) = Y[(X- lnj ), 

3=1 

n . n+1 

B»W = E II - 7^i,r - y) e- 4 ^ , (4.9) 

j— I s j-j j r= i 
n , n— 1 

<y») = - e n n<^ - iw-* + 1 ) ««« 

define a representation of the Yangian. 
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We notice that the explicit expressions for A n (X), B n (X), C n (X) may be obtained directly 
from the defining relations of the Yangian. One may start with the following representation 
of the maximal commutative subalgebra A(qI(N)) by polynomials with real roots: A n (X) = 

n 

n (A — 7nj), n = 1 . . . , N . Then, we should resolve the rest of the Yangian relations to find 
i=i 

the explicit expressions for the generators B n (X) and C n (X) in terms of some operators acting 
on the space of functions depending on the variables j n j ,j = 1, . . . , n; n = 1, . . . ,N— 1 . The 
full set of Yangian relations in terms of A n (A), B n (A), C n (A) is known implicitly through 
the "new realization" by Drinfeld [22j (example after Theorem 1). This may be used to fix 
B n (A),C„(A) up to conjugation by an arbitrary function of j n j . The expressions in (14. 9J) 
are obtained by choosing the appropriate representative. 

The following integral formulas (which may be considered as a "quantization" of the clas- 
sical relations (|3.18|0 express the generators of U(gl(N)) in terms of A n (X) , B n (X) , C n (X) : 



E n , n+1 = A-\X)B n {X)dX 



E n+ x, n = 1 ^jc n {X)A- 1 {X)dX, (4.10) 
1 £ A n (X) dX 1 



E nn = (b . fe (n — 1) . 

2uhJ X-i(A-f) A 2 l ; 

Here, the integrands are understood as Laurent series and the contours of integrations are 
taken around oo . This representation is similar to the expressions of the generators in the 
"new realization" of the Yangian through A n (A), B n (A), C n (A) [H2] (example after Theorem 



5 Application to quantum integrable systems 
5.1 The open Toda chain 

The open Toda chain corresponding to gl(N, R)) is one of the simplest examples of an 
integrable quantum mechanical system with iV degrees of freedom (see [10] and references 
therein). It has iV mutually commuting Hamiltonians, the first two of which are given by 

TV 

i=i 

(5.1) 

N-l 

j<k j=l 

where [x n ,p m ] = ih5 nm . 

Recall the general idea behind the Representation Theory approach to quantum integra- 
bility in the case of the open Toda chain. Let V and V be any dual irreducible Whittaker 
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U(qI(N, R))- modules and id6^ and w' G V be the corresponding cyclic Whittaker vec- 
tors with respective characters x+{En,n+i) — —ihr 1 and X-(En+i,n) = —iff 1 . We assume 
that the action of the Cartan subalgebra is integrated to the action of the Cartan torus, so 
that the following function is well defined 

JV 

/ -tf-pW/ / -p** E n v (5.2) 

This is nothing but the GL(N,M) Whittaker function [S] written in terms of the Gauss 
decomposition ^Tj. In the case of irreducible Whittaker modules, the action of the elements 
of the centre Z of the universal enveloping algebra U(qI(N, R)) is proportional to the action 
of the unit operator [S] (Theorem 3.6.1). Consider the function 

N 

T _ -x-pW/ i (5.3) 

where z belongs to the centre Z . Using the properties of the Whittaker vectors one could 
show that there is a differential operator V z in the variables Xk such that 4 , lN1 ,...,' 7NN = 
T> z ip JNlr „ )1NN . Thus, taking the first two Casimir operators, one gets the first two Hamiltoni- 
ans (|5.1jl . On the other hand, z is an element of the centre Z and ipf m ,...,y NN is proportional 
.,-jnn with some numerical coefficient. Hence, we get the common eigenfunction for 
the set of differential operators corresponding to the elements of Z . 

Let us give an explicit realization of the representation of U(qI(N, R)) in terms of some 
difference/differential operators, which integrates to the representation of the Cartan group. 
This will lead to the integral formula for the wave function of the open Toda chain. To find 
such a formula for the GZ-representation, we substitute the expressions (|2.9|) . (|2.10|l . and 
(|2.1aj) into (|5.2|) . thus obtaining 



N(N-l) 



n=l 11 I 1 I ifi J| n=l 



2 — *<p i<« 

where by definition 7 nj - = for j > n . 

In the study of the analytic properties of this solution with respect to 7^ , the following 
reformulation is very useful: 

Theorem 5.1 An analytical continuation of the eigenfunction \5.4\} as a function of *y N 
can be expressed in the form 



I N 

n n+1 



S 



N-i n n h—^- t{^^) 4 g ^ (5 . 5) 

n fc=lm=l R ni^l TT/I 

n=l 11 V / n =l 

s¥=P j<n 
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where the domain of integration S is defined by the conditions min{Im7fcj} > max{Ini7fc +1 m } 

3 m 

for all k = 1, . . . , N — 1 . The integral \5.5\) converges absolutely. 

Proof. Let us change the variables of integration in (|5.4|) : 

'ft n 

ln 3 - In) ~ - J2p { s N) ' (n = l,...,N-l), (5.6) 

s=l 

and recall that pi N ^ = |(iV — 2s + 1) . After elementary calculations, the integral in (|5.4j) 
acquires the form (|5.5|) . It is worth mentioning that after the change of variables ()5.6|) 

N(N—i) 

we shift the domain of integration K 2 to the complex plane in such a way that the 
domain of integration over the variables 7n-i,j lies above the domain of integration over the 
variables j n j . Thus, we arrive at the analytic continuation of the wave function described 
in the theorem. ■ 



Example 5.1 Let N = 2 . In this case the general formula \5.5\) acquires the form 

"0721,722 (-^l) ^2) = ^ 



e £( 721 +722)* 2 



•ia+oo 



X 



c?7 



(5.7) 



11 ) 



where a > max{Ini7 21 , Im722} • Using the formulas 7.2(15) and 7.12(34) from 14 1)1 , one 
obtains the solution for the GL(2, R) Whittaker function in terms of the Macdonald function 
K v {z) : 

^ yn ^(x^,x 2 ) = 4irheM^+y™^+^K ^-^ (jj e (**-»*)/ 2 ) . (5.8) 

The explicit form of the eigenfunction and the absolute convergence of the integral lead 
to a recursive relation between the wave functions corresponding to 0t(n) and g[(n — 1) ; 
this relation is a direct consequence of the Gelfand-Zetlin inductive procedure: 

Corollary 5.1 Let us fix the solution -ip 7ll (xi) = e^ lllXl of the GL(1,M) Toda chain. For 
any n = 2, . . . , N there is a recursive relation between the wave functions ip~ in and : 



1p~f n \ x li ■ ■ ■ 1 x n) 



^ J (7„-l)/n(7„,7„-l)e i=1 fc = 1 ^n-i^l' 



n-1 



Sn 



n-l,j 

3=1 



(5.9) 



where 



^- i) (7„- 1 )=n[ r (^ 



7n— l,i 



1 h 



-1 



(5.10) 
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and 

n—l n 

/n(7n>7„-i) = 11 11 K lh F[ ^ )■ (5-11) 

fc=l m=l 

The domain of integration S n is defined by the conditions min{Im7 n _ij} > max{Im7 nm }. 

j " m 

The above results are in complete agreement with the recursive relations between the wave 
functions and the explicit solution in terms of the Mellin-Barnes representation obtained 
previously in [22] by the QISM-based approach. 



5.2 The hyperbolic Sutherland model 

In this section we give a short derivation of the integral representation for a wave function 
of the hyperbolic Sutherland model. 

The quantum hyperbolic Sutherland model corresponding to Ql{N) is a integrable quan- 
tum mechanical system with iV degrees of freedom (see [2| and references therein). It has 
N mutually commuting Hamiltonians, the first two of which are given by 

N 

n=1 (5.12) 
h 2 /A 



h=J2 {PnPm + —TTf 

z — ' i. sinn ix™ — x„) ■> 



The eigenfunction of the Hamiltonians of the quantum hyperbolic Sutherland model is de- 
fined by the equations: 



N 



n=l 

^ < smli (x m - x n ) J 



(5.13) 

h 2 /A 



In a similar fashion to the open Toda chain, the hyperbolic Sutherland model can be 
explicitly solved using the Representation Theory approach. Let V be a representation of 
U(qI(N, R)) and v N be a vector in V such that 

{E n>n+1 -E n+1>n )v N = 0, (n = l,...,N-l). (5.14) 

Consider the matrix element $(xi, . . . , x N ) of the form 

N 

rfi / -£*> E >* v (5.15) 

$ = (v N ,e k ^ v N ) . v > 
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One can show that ()5.15|) gives (up to a simple factor) the eigenfunction of the Hamiltonian of 
the hyperbolic Sutherland model. To derive this, consider the more general matrix elements: 



N 



- E *kE kk (5.16) 
<P r = (v N , e a Nr v N ) , v > 

where the coefficients of the expansion of *4jv(A) (see fJ2.17jl ): A N {X) = \ N — 

X N ~ 1 a N i + X N ~ 2 a N2 + • • • • The first two have the form: 

JV 



a N i = ih s y^ j E nn , 

71=1 

«jv2 = — ft 2 j ^ E mm E nn — 'E nn — ^ E nm E mn + er 2 (p^)}- 



m<n n=l m<n 



This should be compared with the discussion of the open Toda chain (|5.3|) . The same 
line of reasoning shows that the matrix element ()5.15|) satisfies the differential equations 

JV 

n=1 (5.18) 
-h 2 { d -9 Xn ~ \ coth(a; m - x n )(d Xm - d Xn )}$ = (tr 2 ( 7jv ) + ftV 2 (p^))$ . 

m<n 

Actually, $ is the GL(N,M.) zonal spherical function 
Finally, the wave function \1/ defined by 

= $ Yl smh 1/2 (xj - x k ) ^ 19 ) 

j<fc 

satisfies the equations (|5.1Hj) . 

Now, we give a new explicit integral formula for the eigenfunction of the hyperbolic 
Sutherland model using the GZ representation discussed in section 2. First, let us find the 
vector v N satisfying (|5.14j) . 



Lemma 5.1 Let 

n n+1 

0n(7n,7„+i) = 11 11 ^ 2in +4 r( +-)■ (5.20) 



2ih 

k=l m=l 



Then the vector 



N-l n 

iK E (n-l) E 7n, 



JV-1 



*'v=< ' ' J"l0n(7n,7n + l) ( 5 ' 21 ) 



2ft 

V ~~ ' 

n=l 



satisfies the equations \5.14 )- 
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Proof. Straightforward check. 



Now we have the following 

Theorem 5.2 The solution to \5.1?A) admits the integral representation: 

j<k 



x 



n ra+1 

N-i n n 



1 1 2i7j " 



B E (7»J-7n-l,j)!Cn 



nfc=lm=l R 4. vmj ,„_i lj; ~ n y-r 
: a e " J= d 7nj ■ 
TT T-l / 7ns -"fnp \ Z J- J. J 

N(N-l) n = 1 11 I 1 I ifi ) I n=l 

R 2 S<P J<™ 

Proof. The only non-obvious thing is the convergence of the integral. However, the integral 
converges absolutely due to the inequality ()2.23|) . ■ 

Example. Let N = 2. Then 

^721,722 ( x i ^2) = sinh 1/2 (2;i - x 2 ) 
x e t(72i+722)x 2 j J^J | r ^ 7n-72m + 1 j| 2 e ^ ll{xi - X2)d ^ ( 5 -23) 

Assuming X\ > x 2 , one can close the contour in the upper half-plane thus calculating the 
integral by residues. The answer is given in terms of hypergeometric functions as follows: 

% 21> , 22 ( Xl ,x 2 ) = sinh 1 / 2 ^ - x 2 ) - ^~ 2h e i(72 1 +722)( a:i+ a ;2 ) 

cosh ^(721 - 722) 

f p/ 721—722 \ 

x e -|(ai-a?a) J c oV(72i -722)^1-^2) v 2ifi / p/1 1 , 721-722 ^ I 721-722 . -2{xi-X2)\ ('5 24") 

I T~Yi 721—722 \ ^2' 2 2ih ' 2i/i ' ) \ ' ) 

< L ^2 lift, I 

p/ 721 —722 \ 

I r -jr(72i -722)^1-^2) ^ 2ift / p/1 1 721-722 1 721-722 . -2(x\-X2Y 
p/1 I 721-722 \ V2' 2 2i/i ' 2i/i ' C 

1 V2 2ifr / 

Up to a trivial factor this solution coincides with the Legendre function P721-722 1 (coshfxi— 
x 2 )) (see jH], formulas 3.2(9), 3.2(27)), namely, 

^721,722 ( x i i^) = sinh 1/2 (xi - x 2 ) 

x - ^ K re i(72i+722)(xi+x2)p TO _^ .(cosh^-xa)). ! 

cosh 2^(721 - 722) 2 ' h 2 

Note, in the limit xi — x 2 00 the wave function (J5.23|) has the asymptotic behaviour 

% 21>722 ( Xl ,x 2 ) ~ c( 7 2i,722)e^ {72ia;i+ ^ :E2) + c( 7 22,72i)e* (722a;i+72ia;2) , (5.26) 
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where 



c(Ai,A 2 ) 




Ai-A 2 \ 
2ih > 



(5.27) 



2ih 



is the Harish-Chandra function. 

Remark 5.1 Since the integrand in is written in the factorized form, it is clear 

that there exists a recursive relation, similar to i5.y\) . between the n and n — 1 particle 
eig en functions of the Sutherland model. This fact is the direct consequence of the Gelfand- 
Zetlin inductive procedure. 

6 Connection with the QISM-based approach 

In section 5 we have outlined the derivation of the integral representation of the wave function 
of an open Toda chain based on the GZ representation of U(gi(N, R)) . As was stressed in 
section 4, this representation bears a direct relationship to Y(gl(N, R)) . This provides us 
with another point of view on the way to obtain the wave functions for the open Toda chain. 
In particular, the algebra generated by the images of the polynomials A n (X) , B n (X) , C n (X) 
under the homomorphism n N : Y(gl(N)) — > U(gl(N)) can successfully be used to derive the 
recursive relations for the wave functions. 

We begin with a short summary of the QISM-based approach to the GL(N, R) Toda 
chain due to Sklyanin |20j. The model is described with the help of the auxiliary 2x2 
matrix Lax operator with spectral parameter: 




(6.1) 



where p n = —ihd Xn . 

The operator (JbM)) satisfies the usual quadratic relation 

R 2 (X - fJ,)(L n (X) ® /)(/ ® LM) = (1® LM)(L n (X) ® I)R 2 (X - (i) 

with the rational 4x4 R -matrix (see ()4.2j0 . The monodromy matrix 



(6.2) 



/ A n (X) B n (X) \ 

(X)=L n (X)...L 1 (X) = 

V C n A D n (X) 



(6.3) 



also satisfies the relation: 



R 2 (X - n)(T n (X) ® /)(/ ® T n (/i)) = (J ® TM)(T n (X) ® I)R 2 (X - //). 



(6.4) 



In particular, the following commutation relations hold: 



[A n (X),AM] = [B n (X),B n (ji)\ = [C n (X),C n (/i)}=0, 



(A- n + ih)A n (X)B n ( f i) = (X-fi)BMA n (X)+ihA n (n)B n (X), 
(A - fi + ih)A n (ii)C n (X) = (A - fi)C n (X)An(fi) + ihA n (X)C n (/j) . 



(6.5) 
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The polynomial A n (X) is a generating function for the quantum Hamiltonians of the GL(n, M) 
Toda chain, and the wave functions ip~ fn (xi, . . . , x n ), n = 1, . . . , N are defined as solutions 
of the differential equations 

n 

A*(A)V> 7nl ,..., 7nB = IJ( A ~ ^j) V>7»l,->7»n ■ (6.6) 
3=1 

The idea of the QISM approach is to solve iteratively these equations for ip^ , n = 1, . . . , N 
using the additional operators B n (X) and C n (X) . From the results of [22] (eq. (2.7)), the 
following expressions for the operators A n (X), B n (X), and C n (X) can be obtained: 

n 

A n (X)=l[(X- lnj ), 

3=1 
n , _ 

W=. I+ "EIlfv^-. (6-7) 

3=1 s+] ■> 

C n (A) = z 1 -" V FT A 7ws e*"** . 

■ , , . ini Ins 
3=1 s^j J 

A simple check shows that, thus defined, the operators (|6.7|) satisfy the relations ()6.5|) . Hence, 
one establishes the connection between the generators ^4„(A), <6 n (A) and C n (X) entering the 
description of the GZ representation (|4.9|) . and the corresponding QISM generators ([6.7)1 : 

A n (A) = A n (X) , 
5 n (A) = s' 1 oi3 n (A)o Sn , (6.8) 
C n (X) = s~\ o C n (X) o s n _ 1 , 

where s n is defined by (12. 10(1 . Therefore, the expressions ()6.5|) should be compared with 
(S3J). 

At the end of this section, we illustrate the connection between the differential operator 
A n (X) entering in ()6.3|) and the generator *4 n (A) defined by (j2.17|) . First, note that (j6.1j) . 
(16.3)1 define the relations between the operators A n (X) for different levels: 

A n (X) = (X-pnjA^tiX) - e x ^- x "A n . 2 (X) , (6.9) 

where n — 1, . . . , N and A_i — 0, A — 1 . 

Proposition 6.1 Lei if)y N1 ,...,y NN be defined by Then the operators A n {X) and A n (X) 
are related by 

N 

A f\\ I -x-pW / I -p x * E kk , i(N — n)h\ v (6.10) 

K{X)v lm ,..., lNN = e p K,e fc=1 X(A- 2 )w JV ). v 1 
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Proof. We prove the statement by showing that the operators A n (X) defined by (|6.10j) 
satisfy the recursive relations (|6.9|) . Consider the relation (|2.19|) . The generators E nn and 
E n)n _ k+ i commute with A n -i(X) and A n -k(\) , respectively Therefore, due to (|2.27|) and 

(£3), 

JV 

(w' N ,e~^ XkEkk A n (X-^i^)w N ) 

N 

= (A - ihpW + ihd Xn ){w' m e ^ XkEkk A-i(A - ^±m )wN) (6-H) 



^-^(w' N ,e'^ XkEkk A n ^(\ - M^m )wN ). 



We leave a more systematic discussion of the interrelation between the two sets of oper- 
ators A n (\) , £> n (A) , C„(A) and A n (\) , B n (\) , C n {\) for a future publication. 
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